Magnetic ordering in GaAlAs:Mn double well structure 
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The magnetic order in the diluted magnetic semiconductor barrier of double AlAs/GaAs: Mn 
quantum well structures is investigated by Monte Carlo simulations. A confinement adapted RKKY 
mechanism is implemented for indirect exchange between Mn ions mediated by holes. It is shown 
that, depending on the barrier width and the hole concentration a ferromagnetic or a spin-glass 
order can be established. 
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In Gai-zMn^As jjj-fjj, a new prototype of Diluted Magnetic Semiconductors (DMS) (|||, the Mn 2+ cations have 
the 3d shell partially filled with five electrons, in such a way that they carry a magnetic moment with S = 5/2. Besides, 
the Mn ion binds a hole to satisfy charge neutrality, what is, in itself, a complicate impurity problem. Two Mn ions 
occupying nearest neighbors positions interact via an anti-ferromagnetic coupling of their magnetic moments. In the 
fee alloys, these interactions are known to be frustrated (see, e.g., Ref. [||), establishing the possibility of settling 
a spin glass phase at low temperature. However, an indirect Ruderman-Kittel-Kasuya-Yosida (RKKY) exchange 
provided by the sp-d interaction between the Mn 2+ spins and the spins of the Fermi gas (unbound holes, or holes 
in an impurity band) competes with the nearest neighbors anti-ferromagnetic interaction. At low magnetic ion 
concentration the indirect exchange mechanism may become the dominant interaction, leading either to a spin-glass 
phase or a ferromagnetic order. 

Recently some groups |l]-|7]] succeeded in producing homogeneous samples of Gai_ x Mn x As alloys with x up to 
7% avoiding the formation of MnAs clusters by using low temperature (200 — 300° C) MBE techniques. Besides 
its practical importance, this kind of DMS introduces an interesting problem from the physical point of view: Mn 
in the alloy is a strong p dopant, the free hole concentration reaching even p = 10 20 ~ 21 cm~ 3 jjj-fj]]. At small Mn 
concentrations, the alloy is a paramagnet and an insulator. As x increases it becomes ferromagnetic, going through 
a non- metal/metal transition for higher concentrations (x ss 0.03), and keeping its ferromagnetic phase. For x above 
7%, the alloy becomes a ferromagnetic insulator. In the metallic phase, the ferromagnetic transition is observed in 
the range of 30 — 100K, depending on the value of x. The ferromagnetic order in the metallic phase is understood as 
resulting from the indirect exchange between the Mn ions mediated by the hole gas. In quantum wells, it seems to 
exist a threshold for the width of the magnetic layer in order for a ferromagnetic phase to appear |l(| . 

In this work we perform Monte Carlo calculations to study the magnetic order resulting from the indirect exchange 
between magnetic moments in a particular symmetric double quantum well DMS structure formed by two GaAs wells 
of width L separated by a Gao.65Alo.35As : Mn DMS barrier width d. In that DMS the Mn 2+ ion concentration 
is taken as x = 0.05. The Mn 2+ ions substitute the cations elements, each of them providing a hole. In addition, 
the DMS is assumed to be in a metallic phase, but the density of free carriers (holes) is only a fraction r of the 
magnetic ion concentration, what is in agreement with experimental data of Ref. |J. A confinement-adapted RKKY 
jnj formalism is used to obtain the indirect exchange for the double quantum well structure: 
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with 4> n (z) representing the eigenfunctions of the potential well, / the sp-d interaction p2|-^4|, and A the normalization 
area for the otherwise free motion in the (x, y) plane. The coordinates (Ri, Zi) describe the position of the impurity i 
in the plane (x, y), and in the growth direction inside the barrier. ^™>™ ((f) is the equivalent to the Lindhard function: 
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It is worthwhile to mention that the hole system is expected to show a spin polarization in a DMS magnetically 
ordered phase. This effect, which is important in spin resonant tunneling experiments, does not result into a major 
change in the order of the Mn ++ impurities. 

The intra-subband contribution of an occupied subband n to the exchange reads: 
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where m* t is the transversal effective mass, and k^ is the n-th subband Fermi wave vector. 
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The contribution of the inter-subband terms cannot be expressed easily in a closed form. Starting over from Eq. 
a) we arrive to: 
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where we used 
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and A„/ jn = 
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E n )/h 2 . The integral in Eq. (^) is, then, straightforward: 
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It is well known [ pd|jl5[ | that, in the mean field approximation, a DMS quantum well with infinite barriers has 
no inter-subband contribution to the Curie- Weiss temperature. For finite barriers and low carrier concentrations, 
these contributions are also small. That is not the case, however, when a DMS layer is inserted in the middle of a 
non-magnetic quantum well jnj, as in the present structure. 

A Monte Carlo (MC) simulation was performed to investigate the magnetic order in the DMS layer. The spin 
sites in that layer belong to the cation fee sublattice. They are distributed randomly and are occupied by Mn 2+ 
ions, with a concentration x. The calculation is performed in a finite box whose axes are parallel to [100] directions, 
of dimensions L x — L y , and L z = N a/2, where a is the lattice parameter of GaAs, and N the number of DMS 
monolayers (ML) in the barrier. Periodic boundary conditions are imposed in the (x, y) plane. Lateral dimensions 
are adjusted in such a way that the total number N s of spins is about 4400, for all L z . Their initial orientations are 
randomly assigned. The energy of the system due to RKKY interaction described by the above 's is calculated, 
and the equilibrium state for a given temperature is sought by changing the individual spin orientation according to 
the Metropolis algorithm jl6) . A slow cooling stepwise process is accomplished making sure that thermal equilibrium 
is reached at every temperature. The resulting spin configuration is taken as the starting configuration for the next 
step with a lower temperature. For every temperature, the average magnetization < M > and the Edwards- Anderson 
order parameter q are calculated. The latter is defined as 



where the expression under the square root is the MC time average. 

We present results for the eight samples described in Table 1. L was chosen typically 6 nm. The DMS barrier width 
was varied from 4 ML to 18 ML and the density fraction r was taken equal to 0.1 and 0.25. In Fig. 1 the normalized 
magnetization < M > is plotted versus temperature. Samples # 3, # 5 and # 8 show a ferromagnetic order at 
T c «35, 55 and 15 K, respectively. The complete saturation at T = K is not achieved, presumably due to boundary 
effects in the z-direction. It is striking that the magnetization curves in those samples are far from the canonical 
Brillouin function, as already experimentally remarked in Ref. OJ , being instead very similar to the ones obtained by 
these authors. On the other hand, samples # 2 and # 7 are still paramagnetic at T « 1 K. The Edwards- Anderson 
order parameter q given by Eq. || is shown in Fig. 2 as a function of temperature. The ferromagnetic samples have 
nearly the same q{T) as < M > (T), what is not surprising. However, samples # 1, # 4 and # 6 show large g-values 
at low temperature, indicating a spin-glass like magnetic order. In particular, samples # 4 and # 6 show a non 
negligible magnetization at low temperatures. This points to an occurrence of a canted spin phase. 

The results seem to indicate that two relevant parameters compete in establishing the magnetic order in the DMS 
layer. On one hand, a ferromagnetic order is expected to settle as the layer width is increased, in accordance with 
what observed in Ref. Jul. On the other hand, depending on the hole density, and because of the oscillatory nature 
of RKKY interaction, in addition to the ferromagnetic couplings, antiferromagnetic couplings can be switched on. 
This is the ingredient which together with disorder are the origin of the spin-glass phase. This is illustrated in Fig. 3 
where the RKKY interactions are plotted versus the in-plane distance Rij for samples # 5, # 6 and # 8. One 
sees that for ferromagnetic samples (# 5 and # 8) Jy- is essentially positive, while for sample # 6, a non negligible 
negative part is present, especially when considering the occurrence of pair couplings with large iiy's. Therefore it is 
not surprising that a (canted) spin-glass phase sets in. 

In conclusion, we have shown that depending on the DMS barrier width of a double QW structure and the hole 
concentration, different magnetic phases can be obtained. For applications, favorable spin configurations can thus be 
designed. In particular, the possibility of having a ferromagnetic order in the DMS barrier of a double QW structure 
is important for spin tunneling and resonant spin tunneling in nanostructures P,fL7|. 
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TABLE I. Samples characteristics: The width of each one of the two GaAs well in the structure is L = 6nm; N is the number 
of DMS barrier monolayers (ML), r is the ratio of hole density to Mn density; T C (K) is the magnetic transition temperature 
for F: ferromagnetic, P: paramagnetic, SG: spin-glass phases. The sign (*) indicates a possible canted phase. The calculation 
is performed with the value of Nq/3 = —1.2 eV according to ref. |13|. 
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0.25 


F 


15 



5 



FIG. 1. Normalized magnetization vs temperature for samples indicated in Table 1. 



FIG. 2. Edwards- Anderson order parameter q vs temperature for samples indicated in Table 1. 



FIG. 3. RKKY exchange interaction Jij vs Rij for samples # 5, # 6 and # 8. 
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